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Abstract 

We prove the existence and uniqueness of an equilibrium state with unit mass to the 
dissipative linear Boltzmann equation with hard-spheres collision kernel describing 
inelastic interactions of a gas particles with a fixed background. The equilibrium state 
is a universal Maxwellian distribution function with the same velocity as field particles 
and with a non-zero temperature lower than the background one, which depends on 
the details of the binary collision. Thanks to the iJ-theorem we then prove strong 
convergence of the solution to the Boltzmann equation towards the equilibrium. 
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1 Introduction 

The present paper deals with the linear dissipative Boltzmann equation for hard spheres in- 
teractions, and solves some questions left open in the very recent work by Toscani and Spiga |26|. 
devoted to the investigation of a linear dissipative Boltzmann equation for Maxwell molecules. 

In the last years, the study of kinetic models for granular flows received a significant interest. 
The largest part of this work deals with kinetic nonlinear models based upon generalizations of 
the Boltzmann-Enskog equation. We refer the interested reader to the review articles |11II18|[T^ . 
However, most of the studies refer to inelastic Maxwell particles, both for the driven case or for 
the free case OlSlin]. Such (pseudo-)Maxwellian models enjoy nice mathematical simplifications 
and lead to exact analytical results |14[ [TE\ (see also the recent developments on the inelastic Kac 
model |24[ \2b\ ). At present, only a few papers consider real interactions among grains, which are 
well described by the inelastic hard-spheres models [Tj El ■ 

Despite their importance for practical applications, linear equations for dissipative models have 
been much less addressed. To our knowledge the only progresses on the matter are those of the 
afore-mentioned paper [211 and of R. Pettersson "23" . Linear models describe the time evolution of 
the distribution function /(cc, v, t) of particles of masses m (representing the granular gas) colliding 
inelastically with particles with masses mi of a fixed background. Throughout this paper, the 
subscript (1) will be addressed to the fixed field particles whose distribution function is known 
and is assumed to be a normalized Maxwellian Mi with given mass velocity and temperature. 
Note that, the grains being cohesionless, long-range interactions of any kind are irrelevant. Thus, 
the only model with real physical interest is the hard-spheres model. As introduced in |2fij . the 
evolution of /(a;, v, t) is given by 



dt 



v,t)+v- V,/(x,v,t) 



1 

2^ 



|q- n| 



1 



/(v,)Afi(wO-/(v)Afi(w) 



dwdn. (1.1) 



Here A denotes the constant mean free path, q is the relative velocity, q = v — w. The velocities 
(vj,,Wj,) are the pre-collisional velocities of the so-called inverse collision, which results in (v,w) 
as post-collisional velocities. In the granular setting, the most important feature of the collision 
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mechanism is its inelastic character which induces that (generally) the total kinetic energy is 
dissipated. The constant parameter < e < 1 is called the restitution coefficient and measures the 
inelasticity of the collisions. Whenever e = 1 we recover the usual linear Boltzmann equation (see 
Section 121 for details). 

One of the main feature of this paper is to prove the existence and uniqueness of the (ho- 
mogeneous) equilibrium state of equation Precisely, we exhibit a (non-trivial) distribution 
function AI (depending on the velocity only) such that 

Q(M) = 

where Q denotes the right-hand side operator of (|l.lf) . When e — I, this question is trivial since 
the conservation of momentum and energy implies 

M(v)Mi(w) = M(v,)Afi(wO 

where M stands for the Maxwellian distribution with same drift velocity and temperature as 
Ml but corresponding to particles of masses m (see Section |21 for more details). Then, one sees 
immediately that the integrand of Q vanishes for / = M. Clearly, for < e < 1 this is no more the 
case and it appears difficult to determine explicitly the (eventual) equilibrium state. Actually, the 
two following questions are far for being trivial: 

• Does an equilibrium state exist? 

• If it does, is it given by some suitable Maxwellian distribution? 

In this paper, we answer positively to both questions showing that the equilibrium state is a 
Maxwellian distribution with the same mean velocity as Mi and with a universal non-zero temper- 
ature lower than the given background temperature. Moreover, this equilibrium state is unique, 
provided its mass is prescribed. This Maxwellian distribution coincide with the one derived in the 
pseudo-Maxwellian case ^21- Actually, it is also possible to show that, as in the non-dissipative 
case, the equilibrium state is universal in the sense that its temperature does not depend on the 
collision kernel, but only on the inelasticity parameter e which determines the details of the binary 
interaction. Let us mention here that our results answer to some open questions from '26' and 
complete the study of 23, where the existence of an equilibrium state was used as an assumption 
for some of the results. 

The two previous questions, as well as the problem of the rate of convergence towards equi- 
librium, have been recently addressed in '26" for the pseudo-Maxwellian approximation. This 
pseudo-Maxwellian approximation consists in replacing the relative velocity q appearing in the 
collision kernel |q ■ n| of Q by the unit vector in the direction of q. The pseudo-Maxwellian model 
enjoys in particular two fundamental properties. First, the associated moment equations are closed 
with respect to the moments of the distribution function. Hence, it is possible to derive the time 
evolution of the drift velocity u{t) and the temperature T{t) of /(v, t) and to predict the mass ve- 
locity and temperature of the eventual equilibrium state. Moreover, as pointed out by A. Bobylev 
[2], Maxwell models lend themselves to a convenient Fourier analysis. These two important prop- 
erties enabled to determine the Maxwellian equilibrium state for the pseudo-Maxwellian model 
and to prove also exponential convergence of the solution to towards the equilibrium (in the 
homogeneous setting). 

Unfortunately, these two tools do not apply to the hard-spheres model Hl.l(l and one has to 
proceed in a different way. The main problem is actually to predict what should be the eventual 
steady-state. To this aim, we derive formally a linear Fokker-Planck equation which is naturally 
associated to the dissipative Boltzmann model Hl.l|l (Section^ through the asymptotic of the graz- 
ing collisions (see e.g. JP", Chapter II.9]). There is a noticeable amount of results on the matter for 
the elastic (nonlinear) Boltzmann equation. We mention here the papers ^22, J£, 22s JJ, which en- 
lighten the connection between the nonlinear Boltzmann equation and the Landau-Fokker-Planck 
equation and the works |16[l21j which describe the same procedure for linear problems. We empha- 
size the fact that the main goal of this analysis is not to prove rigorously any kind of asymptotic 
procedure, even if it is possible to make our results rigorous (see Appendix). The method of Sec- 
tion 3 must only be viewed as a formal (but efficient) way to find a suitable approximation of the 
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collision operator Q which maintains the equilibrium distribution. To our knowledge, it is the first 
time (in kinetic theory) that such a limiting process is performed within this scope. The utility of 
this procedure comes out from the fact that the equilibrium state of the Fokker-Planck operator 
is immediate to obtain. It is a Maxwellian distribution which appears then as the candidate to 
be the stationary solution of The main problem is then to show that this Maxwellian is 

effectively a steady state for Q. This will be done by means of Fourier transform arguments. 

A second important task addressed in this paper is the large-time behavior of the solution to the 
space homogeneous version of (jl.lf) . Actually, once the existence and uniqueness of an equilibrium 
state established, the question of trend towards this equilibrium is of primary importance in kinetic 
theory. In this paper, we are able to show a strong L^-convergence result to the Maxwellian steady- 
state. As in , our result is based upon weak-compactness arguments by means of the so-called 
i?-theorem and some estimates on the second moment of the solution to 1)1.1(1 . 

Let us explain now in some details the organization of the paper. In Section 2, we describe 
briefly the dissipative Boltzmann linear model and its properties. In Section 3 we deal with the 
asymptotics of the grazing collision, from which we derive the Fokker-Planck approximation of Q 
which help us to identify the (possible) equilibrium state. Then in Section 4 we show that the 
Maxwellian obtained through the grazing collisions procedure is really a stationary solution to 
(11.1(1 . Thanks to the Boltzmann TJ-theorem we then prove in Section 5 that the equilibrium state 
is unique (provided its mass is prescribed) and that the solution to 1(1. 1() converges (in the strong 
L^-sense) towards the equilibrium. Finally, we end this paper by presenting some open questions 
and perspectives. 



2 The dissipative linear Boltzmann equation 

As briefly described in the introduction, in this paper we are concerned with the evolution of 
the distribution function /(v,t) of granular gas particles with masses m which undergo inelastic 
collisions with the fleld particles (of masses mi) of a fixed background. The background is supposed 
to be at thermodynamical equilibrium with given temperature Ti and given mass velocity Ui, i. e. 
its distribution function is the normalized Maxwellian 

/ TOi \^^^ , mi(v — Ui)^^ 
^^^(-)=(dTj ---P^- 2Ti > 

The main feature of the binary dissipative collisions is that part of the normal relative velocity is 
lost in the interaction, so that 

(v* - w*) • n = -e(v - v^f) • n, (2.1) 

where n G is the unit vector in the direction of impact, (v,w) stand for the velocities before 
impact whereas (v*, w*) denote the post-coUisional velocities. The so-called (constant) restitution 
coefficient e is such that < e < 1, the case e = 1 corresponding to elastic collisions. Thanks to 
(|^ and assuming the conservation of momentum 



one then finds the collision mechanism 

('v*=v-2a(l-/3)[(v-w)-n]n 

|w* = w + 2(1 - a){l - /5)[(v - w) ■ n]n, 

where a is the mass ratio and (3 denotes the inelasticity parameter 

mi a 1^*^ 



m + mi 



The parameter a is such that < a < 1 (we exclude the limiting cases of Lorentz and Rayleigh 
gases), while the inelasticity parameter satisfies < /3 < 1/2. We refer to ^21 for a detailed 
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description of the geometry of the coUisions. It is easy to see that system (12.211 is invertible and 
provides the pre-colhsional velocities of the so-called inverse collisions, resulting in (v, w) as post- 
collisional velocities 

1 - /3 r/ X 1 

V* = V - 2a ^ _ [(v - w) ■ n]n 

1-/3 

= w + 2(1 - a) ^ _ 2^ [(^ ^ ■ 

In contrast to the elastic case (e = 1), the binary collision dissipates part of the kinetic energy 

mmi 



m|v*p + mi|w*p - (to|vP + mi|wp) = -4- 



-/3(l-/3)|q-n|^ <0. 



m + mi 

In space homogeneous conditions, upon using A as a time scale, equation (|1.1|) can be re-written 
in the dimensionlcss form 



df,, 1 



|q- n| 



1 



/(v,)Mi(w,) -/(v)Afi(w) 



dwdn. 



(2.3) 



The factor in the gain term above appears respectively from the Jacobian of the transformation 
dv^dw^, into dvdw and from the length of the cylinders |q*-n|=e|q-n| (see [TJ for details). Let 
Q be the (dissipative) linear Boltzmann collision operator (acting only on the velocity space) 



[q-nl 



4/(v.)Mi(w.)-/(v)Mi(w) 



dwdn 



(2.4) 



The change of variables n ^ — n leads to the equivalent expression 



Q(/) = -/ i/(q-n/|q|)q.n 



4./(v*)Afi(wO-/(v)Afi(w) 



dwdn. 



where H{-) is the Heavyside step function. We can also define the collision operator by its action 
on the observables. Precisely, for any regular test-function ■(/'(v) 



Ztt 



'xR3xS2 



|q • n|/(v)Afi(w)[?/'(v*) - i/^(v)]dvdwdn. 



(2.5) 



Clearly, '0(v) = 1 is a collision invariant (mass conservation) whereas, in contrast to the elastic 
case, ipiv) — V and ?/'(v) = are not (dispersion of kinetic energy). Note that an important 
feature of the hard-spheres model is that (even in the elastic case) the moments equations of Q(/) 
are not closed with respect to the ones of /. 



3 The Fokker— Planck approximation 

The main goal of this section is the formal derivation of a linear Fokker-Planck equation, 
obtained from (|2.3|) through a kind of grazing collisions asymptotics. We point out that, our aim 
in this paper, is not to prove rigorously the convergence of the (re-scaled) dissipative Boltzmann 
operator Q towards the Fokker-Planck operator Qpp H3.9(l below as the collisions become grazing. 
The limiting process we perform here must only be seen as an efficient tool to predict the nature 
of the equilibrium state of Q (if it exists). Nevertheless, this approximation result can be made 
rigorous and this shall be done in the Appendix. Here we proceed only at a formal level. Let us 
assume that all the collisions concentrate around 

|q-n|/|q|^0. (3.1) 

Consequently, according to H2.2|l one has |v* — v| ~ and, for any smooth function -0, one can 
perform a Taylor expansion of '0(v*) around v leading, at the second order, to: 

^(v*) = V(v) + VvV(v) • (v* - v) + ilD)V(v)(v* - v) ® (v* - v) + o(|v* - vp) 

= ^(v) - 2a(l - /3)(q • n)Vvi^(v) ■ n + (^"(l - ^))q ' ") p2^^(^) . n ® n + o(|q ■ n/\q\\^). 

(3.2) 
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In H3.2|l D^-f/; is the Hessian matrix of ip. The o(|v* — vp) term will be neglected in the sequel. 
One clearly observes that the expansion H3.2|l is similar to that obtained in the study of elastic 
collisions between particles of same masses 0]. This property obviously comes from the fact that 
()3.1fl implies also v' ^ v where v' is the post-collisional velocity in the elastic case. The only 
difference is that, in the classical elastic theory, the multiplicative constant 2a(l — /?) is taken to 
be equal to 1/2. Thus, staying at a formal level, dissipative collisions between particles of unequal 
masses does not lead to supplementary difficulties. 

Let us consider a referential frame with the a;-axis directed along q. Then, 

'n = (cos 6*, sin 6* cos ^, sin 6* sin ^) < 6* < 7r/2, < C < Stt, 

Q nl (3 3) 

cos 9 = — I — I — and dn = sin6'd6'df. 

|q| 

Assuming that the collisions concentrate around 9 7r/2, we define 

bsiS)^lx[./2-s,./2]i0) iS>0) 

and 

i-tt/2 



Is= I bs{9) cos^ 9 sm9d9. 
Jo 



One sees that 

h-TT as 5-0. (3.4) 
Consequently, let us define the associated collision kernel 

Bs{ci,n) = bs{9)\ci-nl 

and denote by the collision operator obtained by replacing |q • n| by 5^'^Bs{c^, n) in H2.3|l . 

Remark 3.1. Note that the introduction of the multiplicative factor 6~'^ can be seen as a suitable 
time-scaling in H2.3|l (see Appendix El for further details). 

Our aim is to show that for small values of 5, the (re-scaled) operator is closed to be a 
Fokker-Planck collision operator. Precisely, let us fix / G L^(M'^) and a smooth test-function 
'0(v). Using Qs into H2.5|l leads to 

(V-, Qsif)) = [ B5(q,n)/(v)Afi(w) [.^(v^) - V^v)] dvdwdn. 

Now, inserting the expansion H3.2|l in the above expression gives the second order approximation 



2a{l-f3) 



/ B5(q,n)(q- n)/(v)Mi(w)VvV'(v) • ndvdwdn 

/ Bsiq,n){q-nf f{v)Mi{-w)B'^tp{v) -iK^ndwd-wdn. (3.5) 

jR^xR^xS^ 



^xR^xS- 

To estimate 3], we first compute the integral with respect to dn. According to (|3.3(l 

rTr/2 



Therefore 



/ B5(q,n)(q • n)ndn = 27r|q| / ^^(q, n)(cos2 6*, 0, 0) sin6' d6' 

J §2 Jq 

I-it/2 

= 27r|qp/ bs{9){cos^ 9,Q,G)sin9A9 
Jo 

= 27r|qp(/5,0,0) = 2^/5|q|q. 
= -Is^^^^^ I |q|V(v)Mi(w)Vv^(v,<) . ^dvdw. 

jR3xR3 |q| 
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Thanks to (|3.4(l . the coefficient in front of the above integral converges to —a{l — /3)/7r as 6 goes 
to 0. Consequently 

|q|2j(^)A^^(^)v^^(^^i). q dvdw. (3.6) 

^ jR3xR3 |q| 



We proceed in the same way for J^. One has first 

/ B5(q,n)(q-n)^n®ndn= |q|M bs{9) cos^ 9 siii0 d0 

rTr/2 



n (g) n d^ 



f ' 1 1 

2n\q\^ / bs{e) cos^ 6 sin 6 Diag[cos2 6, - sin^ 9, - sin^ e]de, 
Jo 2 2 

(3.7) 



where Diag[ai, a2, as] is the diagonal matrix in M.^ x K'^ whose diagonal entries are Ui {i — 1,2,3). 
Now, defining 

f.ir/2 



one gets 



Jo 

[ Bs{q, n) (q • n)' n ndn - 27r|q|3Diag[i^5, his - Kg), his - Kg)], 
JS2 z z 



and 

\2 



j2 = (^"(^ J^)) /■ |q|3/(v)Afi(w)D2v,(v,i) • Diag[A',, 1(7, - Ks), his ~ Ks)]dvd^v 
Since is negligible with respect to Ig one concludes with the following approximation 

j2 ^ "'(^J Pf I |q|3/(v)Afi(w)BV(v) • Diag[0, 1, 1] dvdw. (3.8) 

jR3xR3 



27r 

For any z G K'^ (z ^ 0), let S(z) be the symmetric matrix 

z ® z 



§(z) = Id - 

i.e. §(z) is the projection on the space orthogonal to z. Then, combining H3.5|) . (|3.t)|) and (|3.8|) . 
one obtains for the approximation 

(V', Q5(/)) /■ |q|2/(v)Mi(w)Vv^(v) • -^dvdw 

I" JR3xR3 |q| 

^a^(l-^ / |q|3/(v)Afi(w)BV(v) •S(v-w)dvdw. 

27r jR3xR3 

Now, straightforward computations, using the fact that 2|q|q = Divv(|v — wpS(v — w)), yield 

(V;, QsU)) I clvVv^(v) / |v-w|3§(v-w) {kA/i(w)Vv/(v) + - ^)/(v)VwMi(w)} dw 

ZT^ Jr3 Jr3 

where the parameters k and fi are defined respectively by 

K = a2(l-/3)2; /i = a(l-/3). 

Since the above approximation is valid for any arbitrary smooth function ip, one sees that, as S 
goes to 0, Qs can be approximated (up to the constant l/27r) by the Fokker-Planck operator 

QFp(g)(v) = Vv • / |v - w|^§(v - w) • {KAfi(w)Vv5(v) + (k - ^)5(v)VwMi(w)}dw. 

JR3 



Let us write Qpp in a nicer way. Using the fact that §(v — w) • (v — w) = and 



J 1 



one has 



QFp(.g)(v) = /^Vv ■ / |v - w|3Mi(w)S(v - w) • (Vv.g(v) - _ ui)g(v)l dw 

7m3 I nT, J ^g^^ 



where A(v) denotes the invertible matrix 



v) • { Vvg(v) + -^^^^{v - ui)g(v) 



/ |v-w|^Mi(w)S(v-w)dw. 

JK3 



Remark 3.2. In a different spirit, J. J. Brey et al. |H] derived a hnear Fokker-Planck equation 
from H2.3|l in the hmit of small mass ratio (a 0). We also point out that it is possible to 
consider (see e. g. ^J\) the quasi-elastic approximation of Q assuming that (i <^ 1. We adopt 
here the grazing collisions asymptotic since it preserves the parameters a and (3 (and therefore 
the elasticity) modifying only the geometry of the collisions. Note however that the existence of 
an equilibrium state for the linear quasi-elastic approximation of Q is an open problem to our 
knowledge. 

Now, one recognizes immediately that the above procedure preserves the (eventual) equilibrium 
state. Precisely, if F is an equilibrium state of Q5, then QFp(i^) — 0. This strongly suggests that 
one has to select the candidate for being an equilibrium state of Q as the one of Qff- Of course, 
the interest of the above procedure lies in the fact that this latter is easy to exhibit. Indeed, it is 
obvious from H3.9I) that the unique solution with unit mass to Qpp (5) = is given by a Maxwellian 
distribution with drift velocity Ui and temperature 

r mi (/i - 
m \ kT] 

This suggests the following dichotomy: 

Either Q(/) — has no non-trivial solution. 

Either the unique solution to Q(/) — with unit mass is the Maxwellian: 

A^(v) = (^) exp{- } vGM^ (3.10) 

where 

T* = ipJii^T,. ,3.11) 

At this point, we remark that the above Maxwellian distribution is exactly the equilibrium 
state found in in the study of the pseudo-Maxwellian approximation of (|2.3|l . This supports 
our belief that M is indeed the steady state of Q and that, moreover, it is also a universal stationary 
solution (independent of the collision kernel) as it occurs in the elastic case. This will be proved 
rigorously in the following section. 

4 Is the Maxwellian the equilibrium state ? 

The problem of finding the equilibrium state of the linear Boltzmann equation H2.3|l has now 
been reduced to determine whether Q(M) = or not, where 

.^ _ , ■■- r "i(v-ui)2 
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with T* . ^^^^fc^T,. Surprisingly, apart for the pecuhar case of the ID-model, we have 
not been able to prove by direct computation that 

Q(M)(v) = 0, V e 

Actually, we prove this result through Fourier analysis, i. e. we show that 

Q(Mj(0 = 0, e e M^ 
where Q(Af)(^) denotes the Fourier transform of Q(Af). By H2.5() . it is given by 

Q{M){£) = I |q ■ n|Af (v)Afi(w)[exp{-i^ ■ v*} - exp{-i^ • v}]dvdwdn. 

It is immediate to show that, up to a translation of the referential frame, one can assume 

ui 0. 

For the sake of simplicity, we introduce the parameter 

/ \ 3/2 

/ mnii ^ 



and we recall that /i — a{l — /3). Then, 



\2TiT* 



Q(M)(0 = / |q|Af(v)Ai^i(w)exp{-i^ vjdvdwx 

X / 1^ -nl (exp{2^^i(q•n)(^n)}-l)dn. 

|q| 



Now, the key point of our computations is the identity 



Af (v)A/i(w) = C cxp [/^w^ + (1 - ^)v2] I 

= C exp {-^ W - + (v - mf] } , 

that follows from the equality — 77 = — (/i — 1) . Then, 

T* /zTi 

[ |q| exp - M)q' + (v - m)'] I exp{-i^ • vjdvdwx 

X / |^-n| (exp{2iAi(q-n)(e-n)}-l)dn. 

|q| 

The change of variables (v, w) (v, q) yields 
Q(Mj(0 = C / |q| exp{-;^M(l - M)qndqx 

JR3 ^M-tl 

X / exp{--^^(v-^q)2}exp{-i^-v}dvx 

X / 1^ -nl (exp{2^^i(q•n)(^n)}-l)dn. 

|q| 

Performing first the second integral leads to 

Q(M)(C) = Cexp{-^e'} / |q|exp{-^(l-Ai)q2}exp{-z/iq-adqx 

X / I rr • n| ( exp{2 i ^ (q • n) • n)} - 1 ) dn, 
JS2 |q| 



where we used the fact that the Fourier transform of the Gaussian cxp{ — (v — u)^/20} is equal to 
Cq exp{— m • f — 0/2^^} for any O > and u G M'^ (here u = /iq and O ~ mi/iiTi) where Cq 
is a muhipUcative constant depending only on O. Now, as pointed out first by A. Bobylev i2j, the 
inner integral on the unit sphere is an isotropic function of the vectors ^ and q and is therefore 
equal to 

/ le-n/ieil (exp{2*/.(q.n)(e-n)}-l)dn. 

Consequently, 

QiMm^CeM-^e} I |q|exp{-^(l-^)qndqx 
2mi J^i 21 1 



/ IC • n/ICI I ( exp{-i M (q ■ C^)} - exp{-i/^q • ^}) dn, 

JS2 



where 

Since the last integral on the unit sphere only depends on |^| and ^ • q, and |^+| — |^|, we conclude 
that 

Q(M)(O = for any ^ e R^. 

We proved 

Theorem 4.1. The Maxwellian distribution 

with T'^ = 1^ T^Ti is an equilibrium state for Q. 

1 — a[l — p) 

Remark 4.2 (Universality of the MaxwelUan) . One can easily generalize the above compu- 
tations to show that the Maxwellian (|3.10() is an equilibrium state of any collision operator Qb 
taking the weak form 

(v., Qs(/)) = / i?(q,n)/(v)Mi(w)[7A(v*) - ^(v)]dvdwdn, (4.1) 

for any smooth function ip. In (|4.1|) the collision kernel _B(-, •) is assumed to be given by 

B(q,n) = |qr6(q.n/|q|) (4.2) 

with — 1 < 7 < 1 and 6(-) nonnegative. This shows that, as it happens in the elastic case, the 
equilibrium state of the dissipative (linear) Boltzmann equation is universal in the sense that it 
does not depend on the collision kernel. 

Remark 4.3. Note that for a general collision kernel i3(q, n) it is convenient to use (|4.1|) as a 
definition for instead of its strong form: 

Qsif)^^ I B(q,n){J(q,n)/(vOAfi(wO-/(v)Afi(w)}dvdwdn 

where the factor J depends of _B, a and /3 in a complicated way. 

5 On the trend to equilibrium 

In this section we investigate the large-time behavior of the solution to the linear dissipative 
Boltzmann equation H2.3|l . Precisely, let /o be a given (nonnegative) distribution fmiction and 
consider the Cauchy problem 

^/(v,t) = Q(/)(v,t) VGR3, t>0 ^j.^^^ 

/(v,i-0)-/o(v). 
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Since the above problem is linear (and homogeneous), it is not difhcult to construct a (nonnegative) 
mild solution to H5.1(l by a simple iterative method. Moreover, this solution is unique and preserves 
the mass 

/ /(v,t)dv= / /o(v)dv foranyi>0. 

For further details, we refer to where the problem is studied in high generality, including the 
spatially inhomogeneous equation with suitable boundary conditions. 

A fundamental task in kinetic theory is to determine whether the solution to (|5.1|l converges 
toward the equilibrium state of Q or not. Assuming the existence of a unique equilibrium distri- 
bution, this result has been proved by R. Petterson [22 for collision kernels of the form H4.2|l with 
— 1 < 7 < 1 (corresponding to hard or soft interactions). 

Actually, the proof of [22 requires to prove that the equilibrium state we exhibited is unique. 
This can be done by means of the iJ-theorem. 

The 7J-theorem. Let us recall here that the linear if-theorem for the dissipative Boltzmann 
equation (|2.3|) has been first established by R. Pettersson (^Hl, assuming the existence of an equi- 
librium state for Q. We point out that, in contrast to what happens in the nonlinear setting, the 
existence of such a steady-state is necessary in order to prove the iJ-theorem. 

Once the existence of a steady state has been established, one can state the corresponding 
H-theorem and to prove the uniqueness of the stationary solution as a corollary. 

We give here an elementary formal proof of the 7?-theorem (see Pettersson |23 for a more 
general result for the linear inhomogeneous equation with suitable boundary conditions). Let 
$ : M+ ^ R be a convex C^-function. The associated entropy functional reads 

i7.(/|Af)./^^M(v)a>(|^)dv, (5.2) 

where Af(v) is the Maxwellian (|3.10() . The iJ-theorem asserts that iJ$(-|M) is a Lyapunov func- 
tional for the solution of the linear Boltzmann equation. 

Theorem 5.1 (iJ-theorem). Let $ : ^ M &e a convex -function and let Jq he a distribution 
function with unit mass such that _ff$(/o|A/) < oo. Then, 

^^H^{f{t)\M)<Q {t>0) (5.3) 

where fit) stands for the (unique) solution to (|5.1() . 
Proof : It is clear that 

>.(m|M,./J,v,o*'(^)a. 

and the proof of H5.3|l amounts to show that 

for any distribution function / with unit mass for which the above integral is meaningful. From 
(EH) 

/ Q(/)(v)$'('^)dv = ^ / |q.n|/(v)Mi(w)x 
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and this last integral is also equal to 



|q • n|M(v)Mi(w) 



/(v) /(v*) 



Af(v) M(v* 



M(v*) 



M(v*) \M{v*)J M(v) V^^(v) 



dvdwdn. 



Now, since (^<i>' ( ^ ) ,Q(Af)) =0, 



Q(/)(v)$' 



M 

/(v) 
M(v) 



dv 



1 

2^ 



X 



|q-n|M(v)Afi(w)x 



/(v) /(v*) 



The conclusion follows since, $ being convex, 

$'(a)(6-a) < $(6)-$(a) 



M(v) M(v*) 

(a, 6 e R) 



/(v*) 
M(v*) 



dvdwdn. (5.5) 



and (ci>(l:),Q(M)) 



0. 



Direct consequence of the above result is the following uniqueness result (due to Pettersson 
[221 in a more general setting). 

Corollary 5.2. The Maxwellian distribution M given by (|XTU|) IS the unique stationary solution 
to H2.3|l with unit mass. 

Proof : Let F be another equilibrium state with unit mass. Then, according to H5.5|l with 

(z - 1)2 
<I>(z) = one sees that 



\q■n\M{-v)Ml{^^) 



F(v) F{yr*) \ F(v*) 



M(v) M(v*) J Af(v*) 



-dvdvifdn = 0. 



This implies that 

/ |q- n|Af(v)Afi(w) 



/R3xR3xS2 



r F(v) F(v*) 
\ A/(v) ~ Af(v*) 



dvdwdn 



SxS2 



|q • n|F(v)Afi(w) 



F(v) F{v*) 
Af(v) ^ Af(v*) 



dvdwdn = 



where this last integral is null since Q{F) = 0. Consequently, one gets that 



A/(v) A/(v*) 
and this last identity leads, as in the elastic case, to F 



for any (v, w) £ 
M. 



The evolution of the second moment. To prove that the solution to the Cauchy problem H5.1|l 
converges towards the (unique) equilibrium state, one has to establish some a priori estimates on 
moments. Actually, in contrast to the elastic case and because of the lack of collision invariants, 
it is not trivial to estimate the evolution of the moments of f{v,t). This difficulty is peculiar to 
the hard-spheres model and does not occurs for pseudo-Maxwellian molecules [221 1 since in this 
case the equations for moments are in close form. For long-range interactions forces, R. Pettersson 
succeeded in proving uniform estimates on the higher moments of /(v,t) (see |23l Theorem 4.1]) 
for various types of kernels. Unlikely, his arguments do not apply to the hard-sphere model. Here 
we show that the second moment of the solution to H5.1|l remains uniformly bounded in time. 
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Precisely, let the initial distribution function /o G i^(K'^) have unit mass and let f{t) be the 
solution to the Cauchy problem (|5.1|l . Recall that for any t > 0, /(v,t) has also unit mass. One 
considers the following second moment of f{t): 



m^^ I /(v,t)(v-ui)Mv. 



Note that T{t) is not stricto sensu the temperature of f{t) which is defined by replacing the velocity 
Ui by the drift velocity of f{t). Define also 



F(t)^ |v - w^P/(v,i)Mi(w)dvdw. 



One has 



F{t) = [ (v-Ui)V(v,i) / Afi(w)dw+ / /(v,i)dv / (w - Ui)2A/i(w)dw 

JR3 Jr3 Jr3 Jr3 

-2/ (v - Ui)/(v,<)dv ■ / (w - Ui)A-/i(w)dw. 

JR3 JR3 



The last integral is equal to zero by definition of Ui . Therefore 



F{t) = -T{t) + — Ti. (5.6) 



Now, from 1)2. 5|l . one has 



dm 



m 



[ |q-n|/(v,t)Mi(w){(v* -ui)2 - (v - Ui)^} dvdwdn 



dt Gn 
and 

(v* - ui)2 - (v - ui)2 = 4a2(^ _ ^)2|q . _ 4^(^ - /3)(q • n)(i; - Ui) • n 

= -4a(l - /3)[1 - a(l - /3)] |q • np + 4a(l - /3)(q ■ n)(w - Ui) • n. 

Consequently, 

^^ = -|^m(1-«(1-/3)) / |q.n|3/(v,i)Afi(w)dvdwdn 

2m /■ 

H /i / |q ■ n|(q • n)((w — Ui) • n)/(v, t)A/i(w)dvdwdn. 

Jr3xR3xS2 

Recall that < = a(l ~ /3) < 1. Since 

/ |q • nj'^dn = vrjqp and / |q • n|(q • n)(w — Ui) • ndn = 7r|q|q • (w — Ui), 

one gets 

dT{t) 2m 



2m f 
<-^/i(l-M)/ |q|V(v,i)Mi(w)dvdw 

2m /■ 

+ |qhw-Ui|/(v,t)Mi(w)dvdw. 

Let us first look for a upper bound to the second integral in terms of F{t). Using the fact that 
/ (w — Ui)|w — Ui| Mi(w)dw = 0, we obtain 

JR3 

/ |qp|w - Ui|/(v,t)A-/i(w)dvdw = / (v - Ui)^/(v, t) / |w - Ui | Mi(w)dw 

JR3XS2 Jr3 Jk3 

+ / /(v,<)dv / |w - Uip7\fi(w)dw. 

JR3 JR3 



12 



Thus, 

/ |q|V - Ui|/(v,i)Mi(w)dvdw < CiF(i), (5.7) 

JR3XS2 



where 



Ci=max<^ / w-Ui Mi(w)dw, -j^ \2T,-r i \a \ 

Uk3 |w-Ui|2Mi(w)dw J 



(5.8) 



is a positive (exphcit) constant depending only on Mi. Moreover, Jensen's inequaUty gives 
/ |q|3/(v,i)Mi(w)dvdw> f / |qp/(v,i)Afi(w)dvdwV . 

jR3xR3 \jR3xR3 / 

Now, combining (|5.7|l and (|5.8|l and (|5.6(l one gets the differential inequality 

^'^^ < -2/i(l - y)F{tf/'^ + 2^CiF(i). (5.9) 



dt 

A direct inspection then shows that the solution to (|5.9|l satisfies the bound 

f(t) <max| ^^^'^^^ ,f(0)| i>0. 
Turning back to T{t) one obtains that 

sup/ (v-ui)2/(v,t)dv < C50 (5.10) 
t>o Jr3 

provided / v^/o(v)dv < co. Note that the above bound for T{t) is explicitly computable in 

iR3 

terms of /o, Ci, a and (i. 

Now, mass conservation and the i/-theorem, together with estimate H5.1U|I imply that, provided 

/ (l + v2 + |log/o(v)l)/o(v)dv<X<cx) (5.11) 
at any subsequent time i > 

/ (l + (v-Ui)2 + |log/(v,i)|)/(v,i)dv<ifi <oo, 

JR3 

and this implies the weak-compactness in L^(M'^) of the family {/(v, i)}t>o. Now, following the 
strategy of 23 , one obtains first the weak-convergence, and using translation continuity, the strong 
convergence towards the equilibrium of /(v, t). 

Theorem 5.3. Let fo G i^(R'^) be a distribution function with unit mass satisfying (|5.11|) and let 
/(v,i) be the solution to the Cauchy problem (|5.1|l . Then 

hm ||/(t)-/o||LMK-^) =0- 

t— »oo ^ ' 

Remark 5.4. We conjecture that, as it occurs for the pseudo-Maxwellian approximation the 
decay of \\f{t) — /o||li(R3) towards is exponential (with an explicit rate). 



6 Appendix 

Let us make rigorous the derivation of the Fokker-Planck equation we formally obtained in 
Section 121 This can be done a posteriori using the results of the previous Section. We maintain 
here the notations of Section |3| The only difference is that, hereafter denotes the collision 
operator with kernel 

B5(q,n) = 65(0)|q-n|. 
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whereas in Section 01 we considered a re-scaled kernel (see Remark . Actually, in Section 
we were concerned with an approximation of the Boltzmann collision operator whereas in this 
Appendix, we approximate the solution to the Boltzmann equation as collisions become grazing. 
This difference will appear clearly in the sequel. Let /o be a nonnegative distribution function 
fulfilling estimate ()5.11f) and consider the following Cauchy problem 

^(v,<) = Q5(/,-)(v,t) t>0, veM3 ^g^^ 

/5(v,0) = /o(v). 

As we outlined in Sectional problem Htj.lfl admits a (unique) weak-solution which satisfies 



di / /A-(v,i)9t7A(v,i)dv- / /o(v)^(v,0)dv 

= dt f S^(q,n)/^(v,^)A/l(w)[V;(v^^)-V;(v,^)]dvdwdn, (6.2) 

for any ip £ C'2c([0, -l-oo[xR^). In usual notations, C2 c{[0, -l-oo[xR^) denotes the space of all 
functions ip which are continuously differentiable with compact support in [0, -|-oo[ and twice 
continuously differentiable in R'^. We recall that 

Ig:^ bsie)cos^9 sinOde as 5-0. (6.3) 

Jo 27r 

Let us introduce the time-scaling 

gs{^^,t)^fsi^^,S-U). 
Then, considering a test-function of the form "05 (v, t) — tp{v, S^'^ t) into H6.2(l leads to 



dt gs{^,t)dti^{v,t)dv - /o(v)V(v,0)dv 

Jr3 Jts.3 

/ B5(q,n).g^(v,^)Ml(w)[1/;(v^^)-1/;(v,^)]dvdwdn. (6.4) 

The key point of the approximation procedure is the following 

Proposition 6.1. There exists a nonnegative function g : [0, +oo[^ i"'^(R'^) and a subsequence, 
still denoted {gs)s>o such that gs converges weakly in Ll^^{[0, +oo[, L^(U.^)) towards g as S goes to 
zero. 

We leave the proof to Proposition 16.11 to the end of this Appendix and explain now how to 
derive the Fokker-Plank equation from it. Inserting the expansion (|3.2(l into H6.4|l leads to the 
approximation 

dt 55(v,t)atV(v,t)dv- / /o(v)V'(v,0)dv= / J]{t)dt+ Jj{t)dt + Rs 

Jo JR3 Jo Jo 

--^4^74^/ / B5(q,n)(q-n)(75(v,i)Afi(w)Vv?A(v,t) -ndvdwdn 

^T^f) Jo JR3xR3xS2 
(2a(l-/3))2 roo 



/ dt [ B5(q,n) (q- n)^g5(v,i)Afi(w)I])VXv,i) • ndvdwdn + R5. 

Jo JR3xR3xS2 



In the above expression, = Tisitp) is a remainder term obtained from the Taylor expansion of 
V' l|3.2|l (see [21] for details). One can prove as in the elastic case [201 that 

limR5(V') = (6.5) 

<5^0 
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for any test-function ijj e C\^{[Q, +oo[xM^). As in Section 13 one can show that 

r 3l{t)At = r At [ |q|25,(v,i)Afi(w)Vv^(v,<) • A dvdw. 

Jo Jq jk3xr3 |q| 

Consequently, thanks to Proposition 16. II 

hm / J]{t)dt = / / |qpg(v,t)Mi(w)Vv^/'(v,t) • -^dvdw. (6.6) 

^^oJo Jo Jr3xr3 |q| 

We proceed in the same way for to obtain 

Urn / J5(i)di = — !^ / dt |qp5(v,i)Mi(w)ID)V(v,t)-Diag[0,l,l]dvdw. (6.7) 

Then, combining (|6.4|) . H6.6|) and H6.7|l . one sees that the weak hmit g satisfies 
dt / g{^,t)dMv,t)dv - f /o(v)V'(v,0)dv = -;^ / dt / dv 

JR3 JR3 ^TJ" Jo JR3 

VvV'lv,*)- / |v - w|^S(v - w) {KAfi(w)Vvg(v,t) + (k - ^).g(v,i)VwMi(w)}dw, (6.8) 



where k and fi are defined in Section |31 At this point it is not difficult to recognize in (|6.8|) the 
weak formulation of the Cauchy problem 

|(v,t).i-Qpp(,)(v,t) t>0,veM3 ^^^^ 

^g(v,0) = /o(v) 

where the Fokker-Planck collision operator is given by (|3.9|l . It is well-known that problem H6.9|l 
admits a (unique) nonnegative weak solution g. We proved 

Theorem 6.2. There exists a subsequence, still denoted gg, such that 

gs g weakly in LI^^{[0,+oo[, L'^{R^)) 
s — >o 

where g is a weak solution to the Fokker-Planck equation ()6.9|) with initial datum fo satisfying 
(lOTTl . 

It remains now to prove Proposition 16.11 Clearly, it is enough to prove that the sequence gs 
satisfies the uniform estimate 



<5> 

Now the estimate 



sup / (l + (v-Ui)2 + |log55(v,<)|)g5(v,i)dv<oo. (6.10) 

>0,t>0 JR3 



sup / (1 + |log55(v, i)|)55(v, i)dv < oo 

>o,t>o Jr3 



5>G 

follows from the iJ-theorem applied to fs{t). Note that the iJ-theorem turns to be valid for the 
collision kernel Bs{-) since the equilibrium state is universal (see Remark \4.2\i . Now to prove the 
remaining estimate, we proceed as we did in Section [S] to derive formula H5.10|l . We only sketch 
here the main changes. Let 



Then, 

dTs{t) _ md'^ 

dt ~ "evT 

Since 



Tsit)^"^ [ g5(v,i)(v-ui)2dv = ^ / fsiv,d-H){v-n,rdv. 

■J JR3 O Jr3 

/ Bsiq ■ n)/5(v, S-H)Miiw) {(v* - Ui)^ - (v - Ui)^} dvdwdn. 

jR3xR3xS2 

/ |q-npS5(q-n)dn = 27r/5|q|3, 
/ Bs{ci • n)(q • n)(w — Ui) • ndn < 27r/5|q|^|w — Ui 

JS2 



and lim5_»o 2tt5^Is = 1, we can draw the same conclusion of Sectional 
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7 Conclusions 



We proved existence and uniqueness of a collision equilibrium for the dissipative linear Boltz- 
mann equation with general collision kernel. This equilibrium state is a universal Maxwellian with 
the same mass velocity as the field particles background and with a (non-zero) temperature al- 
ways lower than the one of the background (depending on mass ratio and inelasticity). As early 
noticed in [26) . this follows from the combined effects of momentum and energy exchange with 
fields particles on the one side and, on the other side, of energy dissipation in the binary collisions. 

We point out that the existence of a Maxwellian equilibrium at non-zero temperature is of 
primary importance to reckon the hydrodynamic equations for the considered granular flow. This 
can be done through a Chapman-Enskog procedure (see the conclusions of [26)'). 

In space homogeneous conditions, we showed that the solution to the linear dissipative Boltz- 
mann equation converges towards the equilibrium state as time goes to infinity for any initial datum 
with finite entropy and temperature. Unlikely, our convergence result is based upon compactness 
arguments and we have not been able to determine the decay rate towards the equilibrium. We 
may hope that, as it occurs for the pseudo-Maxwellian approximation |26| . the relaxation to equi- 
librium is exponential. The results of Appendix show that the Fokker-Planck equation (|6.9() is a 
good approximation of H2.3|l when collisions become grazing. It is well-known (see 1 ) that the 
solution to H6.9|l relaxes to M exponentially with an explicit rate related to T'^. This supports us 
in the belief that the same occurs for the dissipative Boltzmann equation H2.3|l . We may infer that 
dissipation-dissipation entropy methods should lead to such a result. Work is in progress in this 
direction. 
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